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Abstract 



■r^lj- . We study bound-state effects on the ti production cross section in the thresh- 

^D I old region at hadron colliders. The bound-state effects are important particularly 

at the LHC where the gluon fusion is the dominant subprocess. Due to the for- 
mation of ti resonances in the J = color-singlet channel of gg — > tt and the 
large width of the top quark, the ti invariant-mass distribution peaks at a few 
^ ' GeV below the ti threshold, and it is significantly enhanced over the naive NLO 

j^ ' prediction until several GeV above the threshold. We present predictions of 

the ti invariant-mass distribution which incorporate both the bound-state effects 
and initial-state radiations up to NLO. The bound-state effects would lead to a 
substantial deformation of top-quark kinematical distributions in the threshold 
region. 
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In the forthcoming CERN Large Hadron Cohider (LHC) experiment, top 
quarks will be produced copiously. A huge top-quark pair production event 
sample is considered as a good template for performing various physics studies as 
well as detector calibrations. For this purpose, it is important that we understand 
physics of the top-quark production and decay accurately. Through dedicated 
studies of the top-quark pair production events at the Fermilab Tevatron, it has 
been recognized that top-quark events can be reconstructed with good accuracy. 
For instance, using a top-quark sample in lepton-|-4-jet mode, the lepton helicity- 
angle distribution in top-quark decays has been measured, which requires a full 
kinematical reconstruction of each event. Up to now, a good agreement with the 
Standard-Model prediction has been observed for this distribution \T\. 

At the LHC, top quarks are produced dominantly via gluon-gluon fusion. 
The gluon distribution function is a rapidly decreasing function of its momentum 
fraction x. Thus, a substantial amount of top quarks are expected to be pair- 
produced close to their production threshold in gluon-gluon fusion at the LHC. 
Therefore, contributions of tt resonances may be important in the threshold 
region. The dominant contribution stems from the J^'-' = ^ {L = S = 0) 
color-singlet ti resonance states. An effective operator for producing one of 
these resonances may be written as 

J^gg-.m^^'"'"'G%G';,^cP{ti). (1) 

In contrast, at the Tevatron, the quark- antiquark fusion is the main process of 
top production. Hence, the tt pairs are predominantly in the color-octet J = 1 
state, and we anticipate that resonance effects are not significant. 

The method for incorporating ii bound-state effects has been developed mainly 
in the studies of ti productions in e~^e~ collisions [21 [3]. Formally, in the limit 
where we neglect the top-quark width, Fj — > 0, bound-state effects can be incor- 
porated by resummation of the threshold singularities {as/P)"", where /3 is the 
velocity of the top quark in the tt cm. frame. It is known that, due to the large 
top-quark width, resonance peaks in the top-quark pair production cross section 
are smeared out, and a broad enhancement of the cross section over the entire 
threshold region results from the bound-state effects. 

On the other hand, a distinct feature of tt productions in hadron collisions 
is that effects of initial-state radiation (ISR) of gluons are significant, especially 
close to the top pair production threshold [U O [6]. They can be incorporated 
by resummation of soft and collinear logarithms, which are given in the form 
(a^ In /?)". So far, up to our knowledge, only two papers [21 [10] studied bound- 
state effects on ti productions at hadron colliders that are relevant for a realistic 
top-quark massO In Ref. [9], the ti invariant-mass distributions are computed 
in the threshold region incorporating the leading-order (LO) bound-state effects 
(Coulomb effects) and ignoring other QCD corrections. In Ref. [TO], in the study 
of the total production cross sections of top quarks, the Coulomb effects above 



* Earlier works [71 [8] considered production of sharp toponium resonances. Their phenomenology 
is rather different from that of today's realistic ti resonances, which have much larger decay widths. 



the tt threshold are examined in LO, in comparison to the next-to-leading order 
(NLO) ISR effects. 

In this paper, we study the tt invariant-mass distributions in the thresh- 
old region, incorporating both bound-state effects and ISR effects at NLO. The 
differences of the present work from the previous ones are as follows. As com- 
pared to Ref. [9], we incorporate the ISR effects and the NLO corrections to the 
bound-state effects. As compared to Ref. [101 , we examine the tt invariant-mass 
distributions; we include contributions from the resonances below the ti thresh- 
old; we include the NLO corrections to the bound-state effects. 

In the computation of the parton-level cross sections, the bound-state effects 
and the ISR effects factorize in the LO of the threshold singularities and in the 
leading- logarithmic (LL) approximation of the soft-|-collinear singularities [lUlG]: 



^iSR{s;i -^ f) = K, 



(c) 



dz a(s' 



zs; I 



f)Ft\z) 



(2) 



Here, d{s';i — > /) is the cross section without ISR (i = qq or gg, f = tt in 
various color and J states); aisR^s'ji — > /) is the cross section including the 
ISR effects. The ISR function F^ depends on the initial-state partons (i) and 
the color of the final-state ti system (c = 1 and 8 for color-singlet and octet 
tt states, respectively). K!>'^ denotes the hard- vertex correction factor, which is 

(c) 

normalized as K^^ = 1 in the leading order. 

In this paper, for simplicity of our analysis, we include the ISR effects up 
to NLO. It is known that the differences between the NLO and next-to- leading 
log (NLL) corrections to the ti production cross sections are small, if we ignore 
bound-state effects |10t [6]. Hence, the NLO approximation may be justified 
for our first analysis of the threshold cross sections. The ISR functions for the 
production of various color states up to NLO are given by [8l |12] 



Ft\z) = 6{1 - z) + 



asim 



vr 



/, 



(c) 



JlF_ 



k 



(c) / fJ-F 



2mt 



S{1 



.(3) 



The 0{as) terms reaql 






= 4C7a 




= 4Ca 




-Ca 


AS) ( , ^'F\ 

^'' V '2m J 


= ACf 




-Ca 



r/ln(l 


-z) 


[[ 1- 
Yln(l 


-z) 


[[ 1- 
( ' ^ 


z 
] 


'fln{l 


I' 

-z) 


[[ 1- 
f 1 ^ 


z 

\ 



1 



l-z 

1 

l-z 



l-z 



In 



In 



In 



fJ-F 
2mt 

2mt 



JlF_ 
2mt 



, (4a) 



(4b) 



l-z. 



(4c) 






MF 

2mi 



fc, 



99 



(c) / MF 



2?TT-/ 



-/5oln 
-3Cf In 



j»F 
2mi 

/"F 



2mtJ' 



(5a) 
(5b) 



with /3o = -yCa — '^iT'q- Here, mt and fip denote the top-quark pole mass and 
the factorization scale, respectively (we take fj-p = "n^t)] cxs is the strong coupling 
constant in the MS scheme; Cp = 4/3 and Ca = 3 are color factors, and we take 
the number of light quark flavors to be Uq = 5. The plus-distribution is defined 
in the standard manner. 

Let us turn to the threshold cross sections without ISR, a{s';i ^ /)• As 
is well known, the 5-wave part of the cross sections are most important in the 
threshold region. Contributions of L > are suppressed at least by 0^ ~ a1 
with respect to the leading S'-wave contribution. For individual cross sections, 
leading S'-wave contributions reside in the following channels: 



^ = 99 ■• f = ti{L = 0, 5 = 0, J = 0, color singlet); 
^ = 99 J f = tt{L = 0, 5 = 0, J = 0, color octet); 
^ = 99; / = tt{L = 0, 5=1, J = 1, color octet). 



(6a) 
(6b) 
(6c) 



There is no color-singlet channel for i = qq at LO. The process gg -^ ti{L = 
0, 5 = J = 1) in the color-singlet channel is forbidden due to the angular mo- 
mentum conservation and Bose statistics (Yang's theorem) , and the same applies 
to the symmetric {d"'^'^) part of the color-octet channel. It is also forbidden in the 
anti-symmetric (/ ) part of the color-octet channel at LO; this is because it is 
a nai've-T odd transition and is hence forbidden at LO due to the time-reversal 
invariance of QCD. We ignore qg subprocess contributions, which are suppressed 
by ttsfi"^ as compared to the LO contribution. 



z)]+ of /i^^ 



and fqq are common to those of 



^ The leading (double-log) terms oc [ln(l — z)/{l 
the pseudo-scalar Higgs boson production via gluon fusion process [13] and the Drell-Yan process ^Mj , 
respectively. 



Following the standard framework developed for threshold cross sections [2,_3j, 
the S'-wave cross section for i ^ f including bound-state effects is given by 

a[s;i^f) = \a(s;i^f)\. x — ^^ -. (7) 

The non-relativistic Green functions are defined by 

mt 



iE + iTt)-\-^ + v!^U^) 



G'-''\x-E) = d^{x), (8) 

where E = \/s' — 2mt is the cm. energy of the ti system measured from the 
threshold; r = \x\; Tt is the decay width of the top quark; VQQj-,(r) is the 
QCD potential between the color-singlet (c = 1) or color-octet (c = 8) static 
quark-antiquark pair. On the other hand, Gq{x]E) is the non-relativistic Green 
function of a free ti pair, which is defined via Eq. ([8]) after setting VQQj-,(r) and 
Tt to zero. 

We use the NLO QCD potential, which reads [15]: 



T4g,(.;,B) = C(^)^ 



1 + 
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{2/3o[ln(^Br)+7£] + 4'^} 



(9) 



with 



CW = -Cf. C'fl = '^-Cf, (10) 

for the MS coupling. Here, je = 0.5772... denotes the Euler constant. The 
QCD potential is renormalization-group invariant, and we evaluate the above 
expression at the Bohr scale of /^b = 20 GeV, and with n^ = 5. 

Eq. d?]) incorporates the QCD bound-state effects between ti up to NLO. It 
also incorporates the top-quark decay-width effects on the cross section a{s'; i — > 
/) up to NLO, provided that the 0{as) correction to the top-quark decay width 
is included in Tf [H |16] . 

We correct the gg — > ti tree-level cross sections in the J = color-singlet 
and octet channels using Eq. d?]), since they include all the leading S'-wave con- 
tributions; we do not include bound-state effects in the J > channels. In the 
qq -^ ti tree-level cross section, there is only color-octet J = 1 channel which 
contains the S'-wave contribution, and hence we correct it by using the octet 
Green function via Eq. ([7]). 

The S-wave Born cross sections behave as o" oc a^P near the threshold (/? <^ 
1). Expanding itisr [Eq. ([2])] in q^, we correctly reproduce the dominant 0{al) 
corrections of the NLO cross sections [4J, namely, a^pin /3, a^/31n/3 and a^f3^ 
terms. Furthermore, by introducing the following hard-vertex correction factors 
in Eq. ([2]), we match o"isr to the NLO cross sections up to the a^/3 term: 

i^f) = i + ^iM/,wfM (12) 

vr \mtj 



with 

Here, /^ij denotes the renormahzation scale in the MS scheme. (We take fiR = 
rrif) The logarithmic part of hi in Eqs. <\13\i are independent of the color of the 
tt system. They cancel the renormalization-scale dependence of the Born cross 
sections which are proportional to ag(^/j). The non- logarithmic part of hf^ are 
extracted from the NLO cross sections for the heavy quarkonium productions; 
hgg (1) from the results of Refs. [3 [12]; hgg (1) and hqg (1) from the results of 
Ref. [12]. Numerically they read 49(1) ^ -3.22, hfj{l) ps -0.92, and 4?(1) ^ 

— 1.61. We note that the term —8/9 in hiJ is missing in Ref. ^^O This is due 
to the fact that non-decoupling effects from heavy-quark loops are incorrectly 
omitted throughout the computations in Ref. [I2]ll 

''qq 

[5], which reads h ^ (1) ~ —1.8. The agreement is fairly good. Similarly, we can 
extract a color-weighted sum of hgj{l) from the NLO gg — > tt cross section [4J as 
jhgg (1) + jhgg (1) ~ 0.8, whlch Is qulte different from the corresponding value 

— 1.58 that follows from Eqs. (|13p. Nevertheless, we confirmed that these two 
rather different values are marginally consistent with each other if we take into 
account the claimed numerical accuracy of Ref. [4J||| We will discuss this issue 
further at the end of the paper. 

(c) (c) q 

In general, non-logarithmic parts of k^ and hi both contribute to the a^P 

(c) (c) 

term of the cross section. How to separate non-logarithmic parts of /c^- and hi is 
scheme dependent. In this paper, we choose a scheme such that non-logarithmic 
part of A; is zero; see Eqs. (I5|). 



We can extract h\g {!) numerically also from the NLO qq — > tt cross section 



* The non-decoupling effect of the top-quark loop in the gluon vacuum polarization in the s-channel 
diagram contributes 2[nt(s = 4to() — Iit{s = 0)] = —8/9, where IIj is the top quark contribution to 
the gluon two-point function. The factor of 2 appears since it is a part of the interference terms. 

^We thank the authors of Ref. [12] for confirming this error in their results. We also thank 
M. Czakon and A. Mitov for providing the analytic expression of the a^/? term of the qq -^ tt cross 
section, prior to its publication ^17J , which was crucial for us to identify the source of the discrepancy 
between Refs. [4] and [12]. 

^A precise statement is as follows. The present discrepancy corresponds to a 7% difference in terms 
of ao, which is the only parameter contributing to the a^(3 term in the fitting function in Ref. >4|. We 
have checked that it is possible to shift the value of ao by 7% without altering the fitting function more 
than 1%, which is the accuracy claimed for the fit, if we adjust the remaining parameters appropriately. 



Finally, convoluting with the parton distribution functions (PDFs) , we obtain 
the hadronic cross section 

cr{s;i^f)= dr -^(r) (Tisr(s = rs; i ^ /) (14) 

= K^-) I'dr^ir) j\za{zTs-i^f)F'>''\z). (15) 

The partonic luminosity is defined by 

'-{t;^f) = ^ h^^l dX2 fa{xi, Hf) fb{x2, fJ'P) S{t - X1X2), (16) 



dr 



{a,b} ■ 



where the summation is over {a, b} = {g, g} for i = gg, and {a, b} = {q, q}, {q, q} 
with q = {u, d, s, c, b} for i = qq. In numerical calculations, we use the CTEQ6M 
(NLO, standard MS scheme) PDF parametrization [18] at the factorization scale 
of lip = fnt- We present our results in the form of the tt invariant-mass distribu- 
tion, which is, in principle, a measurable quantity from the final-states at hadron 
colliders: 



(17) 



Here, rrijf denotes the invariant-mass of tt, and tq = m^^/s. The ISR function 

F^'^ [z) is convoluted with the partonic luminosity but not with a, which is 
evaluated at fixed s' = m^p 

Our formulas for a{s'; i ^ f) rely on the non-relativistic QCD framework [3] 
and are valid only in the threshold region. More specifically, our formulas are 
subject to 0{a^(3'^) corrections that grow with E (part of NNLO corrections). 
As a result, the tt invariant-mass distributions described above are valid only 
in the threshold region f3 <^ 1; they approach the NLO results of Ref. |;l9j in 
the region a^ <C /? ^ 1 but deviate from the NLO results at higher invariant- 
mass P ~ 0(1). Hence, our predictions should be smoothly interpolated to the 
corresponding NLO results in the region as <C /3 ^ 1. Numerically smooth in- 
terpolations may be performed at ttt-jj — 2mi ~ 5-20 GeV|j 

Below we examine the ti invariant-mass distributions numerically. We com- 
pare the invariant-mass distributions which include QCD corrections in four dif- 
ferent ways: 



^ In this regard, we note that the ratio of the Green functions in Eq. ([7]) has the following high 
energy behavior: 
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Fig. 1: tt invariant-mass distributions for gg — * tt in the threshold region, in the 
J = 0, color-singlet channel at the LHC The cross sections calculated in the Born 
approximation (black dotted), with bound-state effects but without ISR effects (green 
dashed), with ISR effects but without bound-state effects (blue dot-dashed), and with 
both bound-state and ISR effects (red solid) are plotted. We take mt = 173 GeV and 
Ti = 1.5 GeV. 



Born: 



NLO: 



Gr-Fnc: 



Gr-Fnc X ISR: 



The distribution at the Born level. On the right-hand side 
of Eq. (J17p . a is replaced by the Born cross section; the ISR 
function F^ (z) is set to (5(1 — z); the hard- vertex factor K^^ 
is set to 1. 

The distribution with ISR effects but without bound-state ef- 
fects. The right-hand side of Eq. ([T71) (except dCi/dr) is re- 
placed by its expansion in Os up to 0(0;^). Note that it in- 
cludes the Os/P part of the Coulomb-gluon exchange effects in 

The distribution with bound-state effects but without ISR ef- 
On the right-hand side of Eq. (jlTh . F^ {z) is set to 



fects. 

5{l-z 

Our full prediction, Eq. 



and K^'^^ to 1. 



In Fig. [H we plot the tt invariant-mass distributions for gg -^ tt in the 
J = 0, color-singlet channel at the LHC. The above four cross sections are shown 
(dotted, dot-dashed, dashed, and solid lines, respectively), for mt = 173 GeV and 
Tt = 1.5 GeV. The effects of Tt are included only through the Green function in 
Eq. ([8]). The cross sections with binding effects exhibit a resonance peak below 
the threshold, rritt = 2mt = 346 GeV. Due to the large width of the top quark, 
multiple resonance peaks are smeared out and only the broad IS" peak remains 
in da/drritt', the feature well-known in the total cross section for e~^e~ — > ti. 

The peak position is at -E = rritt — 2mt = —2.1 GeV, which is very close to 
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Fig. 2: Same as Fig. [T], but for gg -^ tt in the color-octet channel and all J states are 
summed. 



the peak position of the singlet Green function IniG'^-'. In fact, the rapid fall off 
of the partonic luminosity dCi/dr reduces the peak position by only about a few- 
tens MeV. On the other hand, the ISR effects scarcely change the peak position. 
This means that the peak positions of the ti invariant-mass distributions in the 
color-singlet channels are, to a very good approximation, the same in e+e~ and 
hadron collider experiments. 

The ISR effects enhance the invariant-mass distribution obtained with only 
the bound-state effects (Gr-Fnc) almost independently of mjf, where the en- 
hancement factor is about 1.2. We also observe that the enhancement by ISR 
effects is much more pronounced without the bound-state effects than after in- 
clusion of bound-state effects. Namely, the ratio of the approximations NLO and 
Born is much larger than that of Gr-Fnc x ISR and Gr-Fnc. This is because a 
substantial part of the former enhancement comes from the ag/ (3 term included 
in NLO, which is genuinely a part of bound-state effects, whereas this term is 
included in the Gr-Fnc approximation. 

Fig. [2] shows the invariant-mass distribution in the gg channel as in Fig. [H 
but for the color-octet states and when all J contributions are summed. We 
have included bound-state effects only in the J = channel. By comparing Gr- 
Fnc x ISR (solid) and NLO (dot-dashed), one finds that the bound-state effects 
are not significant. More precisely, among the bound-state effects (a^//?)", con- 
tributions from n > 2 are rather small for the color-octet state. A small tail of 
the cross section below the threshold originates from the large width of the top 
quark. The ratio of Gr-Fnc x ISR and Gr-Fnc is almost independent of ra^i and 
about 1.4. 

The striking difference of the bound-state effects in the color-singlet and 
color-octet channels can be understood as follows. The force between t and 
i is attractive in the color-singlet channel, leading to formation of resonance 
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Fig. 3: Same as Fig. [T] {gg -^ tt), but the color and J are summed. 



states below the threshold m^f < 2mt. On the other hand, the interquark force is 
repulsive in the color-octet channel and affects only the continuum states weakly. 

Fig. [3] shows a similar comparison of the tt invariant-mass distributions after 
the sum of the color-singlet and color-octet for all J states are taken, for the gg 
initial state. They are obtained as the sum of the corresponding distributions in 
Fig. [1] (color-singlet, J = 0), those in Fig. [2] (color-octet, all J), and the small 
contributions of J > 2 states in the color-singlet channel, which are not shown 
separately. 

Fig. m is for the qq initial state, where only the color-octet J = 1 channel 
exists. Each distribution is similar to the corresponding one in the color-octet 
gg -^ ti channel, see Fig. [2j Note, however, that the qq channel gets smaller 
enhancement by ISR than the gg color-octet channel does (the ratio of Gr- 
FncxISR and Gr-Fnc is about 1.1 for the qq channel, whereas the ratio is about 
1.4 for the color-octet gg channel.). The difference of the enhancement factors 
originates from the different color factors in the ISR functions, see Eqs. ([HHS]). 

Displayed in Fig. [5] are the invariant-mass distributions for tt pair production 
at the LHC, where we show explicitly the individual contributions from the gg 
color-singlet (dashed), gg color-octet (dot-dashed), and qq (dotted) channels, as 
well as the sum of them (solid). Thick lines include both bound-state and ISR 
effects (Gr-Fnc X ISR), while thin lines represent the NLO cross sections. The 
invariant-mass distribution for the sum of all channels still exhibits the 15 peak 
below the tt threshold, while it gradually approaches the NLO distribution above 
the threshold. The color-singlet (J=0) gg channel dominates the cross section 
below and near the tt threshold, while the color-octet gg channel is dominant 
above the threshold. We expect that the cross section below the thresholcO 
(336 GeV < rritt < 2mt = 346 GeV) is about 6 pb, or 6 x 10^ events with 10 fb'^ 



^ The validity range of our formula is given by Im^i — 2mt\ ^ a^rrit, which is of order 5-10 GeV. 
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Fig. 4: Same as Fig. [3], but in qq channel (color-octet and J 



The total enhancement due to the binding effects (—10 GeV < m^f — 2mt < 
10 GeV) is estimated to be about 8 pb, which corresponds to 1% of the total 
NLO cross section. It will be important to take them into account in detector- 
response calibration and in the precise measurements of the top-quark mass. 

Fig.[6]is the same as Fig.[5]but for the Tevatron, pp collision at y^ = 1.96 TeV. 
Due to the dominance of qq channel which is color-octet in the leading order, 
bound-state effects are less significant in the threshold region. The increase of 
the cross section above the threshold due to bound-state effects, as defined by 
the difference between Gr-Fnc xISR and NLO cross sections at < rriii— 2mt < 
10 GeV, is about 0.01 pb. The ratio to the total cross section in NLO [4] is 
2 X 10~^, which is consistent with the corresponding value of Ref. [lOj. However, 
the tt pair-production cross section below the threshold (—10 GeV < mtt — 2mt < 
0) is estimated to be about 0.07 pb, or 70 events with 1 fb~^. Even though their 
contribution is still a small fraction, a few clean events in this region could 
significantly affect the top-quark mass determination at the Tevatron. 

We may summarize the new aspects of our predictions over the previous works 
as follows. The most important difference of our predictions compared to those of 
Ref. [9] is the inclusion of ISR effects. These effects increase the invariant-mass 
distributions (Gr-Fnc) almost independently of m^f by about 20%, 40%, and 
10% for the gg — > it color-singlet, color-octet, and qq — > it color-octet channels, 
respectively. Furthermore, our predictions are more stable against variations 
of fip, AtB or fiR, since we include the NLO corrections to bound-state effects 
and ISR effects, while only LO bound-state effects are incorporated in Ref. [9j. 
Concerning the estimates of bound-state effects on the ti total cross sections in 
Ref. [To], we find that they are significant underestimates, since only the effects 
above the threshold have been taken into account. 



Summing up, we find that the bound-state effects significantly alter the 
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Fig. 5: tt invariant-mass distributions for the top-quark production in the threshold 
region at LHC. The color-singlet (green dashed) and octet (red dot-dashed) in gg 
channel, color-octet in qq channel (blue dotted), and the sum of them (black solid) 
are plotted. Thick lines include both bound-state and ISR effects, while thin lines 
represent the cross sections with only NLO effects. 



invariant-mass distributions of the tf production close to the threshold at the 
LHC. The effects will be particularly important for the determination of the 
top-quark mass [19], as well as when the top-quark sample is to be used for cal- 
ibration of jet energy scale, etc. at an early stage of the LHC operation. In this 
regard, it should be noted (known from the studies of e"^e~ -^ tt [2D]) that the 
bound-state effects distort the momentum distribution of top quarks strongly 
below the threshold (m-tj < '2mt). On the other hand, Tevatron experiments 
are almost free from the bound-state effects, because of the dominance of the 
color-octet qq annihilation contribution, as can be seen from Fig. [6l 

We consider that our calculations incorporate the most important part of the 
bound-state effects and ISR effects to the ti threshold cross sections at hadron 
colliders. There are, however, many corrections that should be included to make 
the predictions more precise. Among them are resummation of the collinear and 
soft logarithms (beyond NLO), non-factorizable corrections, etc. They are ex- 
pected to modify our results at the level of 10-30%. Furthermore, as already 
stated, presently there remains a considerable disagreement in a color-weighted 
sum of the hard- vertex factors hgg . Hence, we examined the following case sepa- 
rately: we use the value of hgg as in Eq. (J13b ). which is determined from the two 
mutually consistent results [HI [12]; on the other hand, we determine the value of 

tt cross section of Ref. [1] is reproduced, 



h 



gg such that the a'lfi term of the gg 



,(8) 



(8) 



i.e. hgg (1) ^ 2.39. With this value of hgg , the cross section in the gg color-octet 
channel is enhanced by about 10% as compared to the case with hgg {!) ss —0.92. 
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Fig. 6: Same as Fig. [5], but for the cross section at Tevatron. 



As a result, the cross section for 'sum of ah' in Fig. [5] becomes more enhanced 
as rrijf increases, where the enhancement is less than 2% below the threshold, 
and is about 6% at rritt = 360 GeV. By contrast, the cross section for 'sum of 
air in Fig. [6] scarcely changes. We also varied hr = fip between mt/2 and 2mt 
and found that the normalization of the cross sections changes by about 10% 
accordingly. 
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